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Abstract To understand and predict hygrobuckling behav¬ 
ior of orthotropic or isotropic wood-based composite 
panels, the closed form equations were derived using both 
the displacement function with a double sine series and 
the energy method under biaxial compressions with an 
all-clamped-edge condition. The critical moisture content 
depended on Poisson’s ratio (v) and was inversely propor¬ 
tional to 1 + v for isotropic panels. It did not depend on the 
modulus of elasticity (MOE) at all for isotropic panels, but 
it did depend on MOE ratios for orthotropic panels. As 
expected, the critical moisture content of plywood was 
twice as large as the that of hardboard owing to the differ¬ 
ence in linear expansions between the two panels. The ap¬ 
plication of optimum thickness and aspect ratios obtained 
by the derived equations could improve hygrobuckling re¬ 
sistance without other chemical treatments that could re¬ 
duce the linear expansion of wood-based panels. This study 
also indicated that it would be better to increase the aspect 
ratio rather than the thickness ratio (alh) from the view¬ 
point economics. 

Key words Hygrobuckling • Orthotropic • Wood-based 
panel • Aspect/thickness ratio 


Introduction 

Most furniture panel components are grid structures, and 
the wood-based panels are fixed on the core with adhesives. 
Among them, the hollow core type of flush door is one of 
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the most important components, as it is relatively expensive 
compared to the side, bottom, shelf, and top components. 
Although the behavior of a wooden flush door due to mois¬ 
ture content change has been studied, most studies have 
focused on predicting the warping magnitude as the mois¬ 
ture content decreases. 1 

As the availability of quality tropical hardwood for ply¬ 
wood production decreased in recent decades, wood-based 
composite panels such as medium density fiberboard 
(MDF), hardboard, and particleboard are being used for the 
same purpose. One of the problems when using these wood- 
based composite panels is their surface waviness. This 
waviness can be observed easily with the naked eye, particu¬ 
larly in the case of doors with surfaces treated with a high- 
gloss coating. It is believed that the waviness is due to 
hygrobuckling because of hygroexpansion of the surface. 
There are few rigorous studies on the hollow core door, 
and its behaviors are understood qualitatively rather than 
quantitatively. 

For isotropic panels in simply-supported and clamped 
conditions, the buckling or thermal buckling has been 
investigated by many researchers for other materials. 2 ' 4 
For example, buckling solutions were reported for ortho- 
tropic panels in a simply-supported condition under biaxial 
compression. 5 For all-edges-clamped conditions, Dickinson 6 
compared the simple solution with the previous one using 
Rayleigh’s method and beam characteristic function. By 
contrast, less attention has been paid to the wood-based 
panel materials developed in relatively recent decades. 

In recent years, thermal buckling of laminated composite 
thick plates has been studied using a finite element method 
(FEM). 7,8 The authors reported the effects of panel geom¬ 
etries, such as each layer’s assembly, thickness ratio, and 
aspect ratio, on thermal buckling. 

Studies on the plate buckling of wood-based panels had 
been conducted on shearing of the wall when an external 
force was applied mainly with uniaxial compression, 9,10 but 
biaxial stress occurs at surfaces via hygroexpansion; hence 
the relevant equation should be adopted. Spalt and Sutton 11 
studied hygrobuckling of strips subjected to a uniaxial com¬ 
pression bar using the equation for critical strain that can be 
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used without buckling and the magnitude of buckling above 
the point of critical stress, derived by Suchsland. 

The main motivation for this study comes from the need 
to understand the surface waviness of wood-based compos¬ 
ite panels, which can be severe, particularly during summer 
when they are stored in the yard. It is important to know the 
critical moisture content where deflection develops rapidly 
with further hygroexpansion. The hygrobuckling of wood- 
based panels such as hardboard and plywood has been in¬ 
vestigated in only a limited manner. In fact, no standards 
are available for selecting an optimum distance between 
cores during manufacture of the hollow core type of flush 
components. 

This study was conducted to look into the effects of the 
material constants and geometry on hygrobuckling of the 
hollow core wooden door subjected to biaxial compressive 
forces as well as buckling by hygroexpansion. In addition, 
the influence of moisture content changes that can cause 
buckling was studied for doors made of plywood or hard¬ 
board under various conditions. 


Theoiy 


In this study, only biaxial forces were considered as 
shown in Fig. 1 because the shear force due to hygro¬ 
expansion is not used on wood-based panels such as 
plywood and hardboard. The energy method was employed 
to solve Eq. (1). 


All-edges-clamped boundary condition 


The critical assumption during application of the energy 
method is the shape of the displacement. Its accuracy de¬ 
pends on the selection of the displacement function. For the 
all-edges-clamped condition, the displacement function 912 
can be represented by a double series [Eq. (3)], satisfying 
the boundary conditions [Eq. (2)]. 


w = dw/dx = 0 for x = 0, a\ w = dw/dy = 0 
for y = 0 ,b 


( 2 ) 


. 7Dc . jiy ^ ^ . rma . tiny 

w = sin— sin— > \a mn sm-sin—— 

a b a b 


( 3 ) 


The buckling mode shape is the same as that of three 
vibration modes, as shown in Fig. 2. Denoting the work of 
external forces as A T and the strain energy of bending as A V 
gives: 


The governing equation for the orthotropic buckled plate 



See Appendix for explanation of the abbreviations and Fjg, j. Coordinate system and dimensions for the biaxial compression 
symbols. condition. See Appendix for abbreviations 



Fig. 2. Buckling mode shapes 
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For m > 1 and n = 1, Eq. (6) is expressed as: 
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Substituting Eq. (3) into Eq. (4), integration by parts leads 
to Eq. (5). It should be noted that all cross terms are zero 
due to the orthogonal properties of buckling modes. 
Therefore, 


In the case of higher modes, m > 2 and n > 2, Eq. (6) 
changes to Eq. (9). 
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The relation between stress and hygroexpansion can be 
represented as Eq. (10). 
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A critical moisture content (MC cr ) can be found by substi¬ 
tuting Eq. (11) into Eq. (6). 


Equating AT to AU, the critical buckling stress by external 
forces when m = n — 1 is expressed as Eq. (6). 


f a 2 > 


4 7T 2 

,2 ® y , 

v * ; 

cr 

37ft 


3Ai 9“ 1A 2 


1 2 (v 21 Ai t 2Ag) 


( 6 ) 


All edges simply supported boundary condition 

Using the same procedures for the edge-clamped condition, 
the critical buckling stress can be found from Eq. (14) if the 
displacement function 13 is represented as Eq. (13) in the 
case of the simply supported condition. 


w = d 2 w/dx 2 = 0 for x = 0, a; w = d 2 w/dy 2 = 0 
for y = 0, ft 


Equation (6) is the same result as with another solution 2 
assuming D n = v 21 D n + 2D M for isotropic panels. 

For m — 1 and n > 1, Eq. (6) is expressed as: 
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Equation (14) is the same as for the other solution. 3 
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Experiments 

The dry process for hardboard and three-ply lauan plywood 
panels was selected for test specimens because these panels 
have been mainly used for furniture: hollow core type flush 
components (door, side, top, and bottom parts). The adhe¬ 
sive for hardboard and plywood was urea formaldehyde 
resin. The thickness ratio (a/h) was adjusted as 400:3 be¬ 
cause the panel dimensions (a X b X h) were 460 X 460 X 
3 mm. Control specimens for the determination of MOE, 
dimensional, and moisture content changes were cut along 
four edges of the buckling specimen (Table 1). Radiata pine 
laminated veneer lumber (LVL) was used for the core ma¬ 
terial with dimensions of 26 X 30 mm. All the specimens cut 
from the panel were preconditioned at 25°C and 60% rela¬ 
tive humidity (RH), but the core specimens were placed 
indoors and coated with oil-based paint to minimize the 
effect of dimensional change due to the RH change. The 
square panels were symmetrically assembled with polyvinyl 
acetate (PVA) adhesive (Fig. 3) and subjected to the con¬ 
stant humidity chamber where the temperature and RH 



Fig. 3. Specimen geometry and assembly setup for test (units are 
millimeters) 


were 25°C and 90%, respectively. Holes of 8mm diameter 
were bored in the middle of the core to reduce the differ¬ 
ence in the equilibrium rate between sample panels for 
buckling and control specimens. Center deflections of the 
assembled panel were measured at time intervals with a dial 
gauge. 

The analytical solutions for clamped and simply- 
supported conditions were compared with the results of 
ANSYS 5.5, a commercial FEM program that has been 
widely used. The element type and mesh size used were 
quadilateral 5x5 elements with eight nodes. They were 
compared with other solutions for thermal or mechanical 
buckling and experimental results. 


Results and discussion 

Validation of the derived solutions 

The derived solutions in this study were compared with 
results of other studies, as shown in Table 2. There are a 
few thermal buckling studies, particularly in the case of the 
all-edges-clamped condition, but the critical temperatures 
for the isotropic square panel ( a/h = 100, a - 2 X 1(T 3 , 
v = 0.3) were the same exactly. Also, the nondimensional 
buckling load, N x b 2 /n 2 (v 2l D n + 2D m ), for orthotropic panels 
subjected to uniaxial force [DJ(v n D n + 2£> S6 ) = 1.543, D n l 
D 22 = 0.321, N y = 0] were a little smaller than Dickinson’s 


Table 2. Comparison of thermal buckling for isotropic and buckling 
for orthotropic panels for the all-edges-clamped condition 

Loading Buckling Aspect ratio Buckling load 

types mod e(m,n) {alb) - 

Reference Present 

Isotropic 


Thermal 

Orthotropic 

1,1 

1.0 

168.71 7 

168.71 

Mechanical 

Uniaxial 

2,1 

1.0 

20.74 6 

(19.74)“ 

18.34 6 

(17.42)“ 

20.45 

compression 

3,1 

1.5 

17.77 


“Values cited from Dickinson 6 


Table 1. Physical and mechanical properties of hardboard and plywood at 25°C, 60% RH 


Property 

No. of 

specimens 

Hardboard 


Plywood 


II 

A 

// 

A 

Moisture content (%) 

32 

6.96 ± 0.36 

(10.65 ± 0.44) 

10.51 ± 0.28 

(13.55 ± 0.32) 

Density (g/cm 3 ) 

32 

0.93 ± 0.03 


0.63 ± 0.02 


Thickness (mm) 

32 

2.98 ± 0.03 

(3.04 ± 0.03) 

2.81 ± 0.01 

(2.83 ± 0.02) 

MOE (GPa) 

16 

5.06 ± 0.17 

4.89 ± 0.58 

11.86 ± 0.59 

3.53 ± 0.39 



(3.55 ± 0.52) 

(3.65 ± 0.55) 

(11.55 ± 0.52) 

(3.22 ± 0.40) 

Coefficient of LE (X10 3 ) 

16 

0.380 ± 0.040 

0.365 ± 0.047 

0.157 ± 0.020 

0.155 ± 0.020 


Numbers in parentheses are values at equilibrium condition (25°C, 90% RH) ± SD 

II, machine or logitudinal direction; A, perpendicular to machine direction; RH, relative humidity; MOE, modulus of elasticity; LE, linear 
expansion 






Table 3. Assumed material constants of plywood and hardboard 


Panel 

£,(GPa) 

£ 2 (GPa) 

G 12 (GPa) 

v, 2 

V 2 1 

o.(mlml%) 

Plywood 

10.0 

2.0 

0.4 

0.12 

0.024 

0.0002 

Flardboard 

4.0 

4.0 

1.6 

0.25 

“ 

0.0005 



Fig. 4. Comparison of Dickinson’s results with present results for iso¬ 
tropic and orthotropic materials by aspect ratios. The relative errors 
were calculated as the percentage of Dickinson’s solution divided bv 
the present solution 


value at two modes, but it was slightly larger than other 
results. 6 

Furthermore, to compare critical moisture contents ob¬ 
tained from this study with other results for different aspect 
and thickness ratios, material constants for usual hardboard 
and plywood were assumed as shown in Table 3. In this 
study, the relative errors of the results were smaller than 
Dickinson’s value with the exception for some modes even 
though errors depended on material properties and aspect 
ratios, as shown in Fig. 4. The difference between the 
present result and those of Dickinson might be attributed to 
one of the displacement functions assumed. However, the 
result of this study might be closer to the exact solution than 
Dickinson’s solution because the energy method always 
causes a slight overestimate compared to that of the exact 
solutions. 

As shown in Tables 4 and 5, the results showed good 
agreement between the present solutions and the FEM 
solution when the thickness ratio ( alh ) was 100 for both 
boundary conditions. As expected, critical moisture con¬ 
tents for the clamped condition were about three times 
larger than those for the simply-supported condition. 
Moreover, the critical moisture content of plywood 
was twice as large as that for hardboard owing to the 
difference in linear expansions between the two panels. 
For the simply-supported condition, the present solution 
was slightly larger than the counterpart of the FEM. In 
contrast, they were smaller than the FEM for the clamped 
condition except for some cases of plywood. It should be 


Table 4. Comparison of critical moisture content in the present 
solution and the FEM solution for plywood and hardboard for the all- 
edges-clamped condition at a constant aspect ratio (alh = 100) 


Aspect ratio 
(alb) 

Plywood 


Hardboard 


Present 

FEM 

Present 

FEM 

0.5 

1.5476 

1.4943 

0.5176 

0.5368 

1.0 

1.6729 

1.6548 

0.7018 

0,7386 

1.5 

2.3341 

2.3263 

1.2248 

1.2714 

2.0 

3.8354 

3.7497 

2.0704 

2.1182 

2.5 

5.3747 

5.5275 

3.2112 

3.2484 

3.0 

7.3658 

7.4991 

4.6321 

4,6433 


Results (critical moisture content) are percents 
FEM, finite element method 


Table 5. Comparison of critical moisture content for the present solu¬ 
tion and the FEM solution for plywood and hardboard for the all edges 


simply-supported condition at 

a constant aspect ratio (alh = 

.100) 

Aspect ratio 
(alb) 

Plywood 


Hardboard 


Present 

FEM 

Present 

FEM 

0.5 

0.4066 

0.3984 

0.1625 

0.1617 

1.0 

0.4643 

0.4570 

0.2632 

0.2565 

1.5 

0.6683 

0.6575 

0.4277 

0.4178 

2.0 

1.0819 

1.0623 

0.6580 

0.6453 

2.5 

1.7273“ 

1.6924* 

0.9541 

0,9384 

3.0 

2.6037* 

2.5448* 

1.3159 

1.2965 


Results (critical moisture content) are percents 
"Buckling occurred at mode (2,1) 


recognized that most of the buckling occurred at mode 
(1,1), but it was not consistent for orthotropic panels such as 
plywood as the degrees of orthotropy and aspect ratio 
increased. This phenomena was compatible with the results 
for the plywood with aspect ratios of both 2.5 and 3.0 for 
the all-edges-clamped condition at which buckling occurred 
at a mode (2,1) in both cases of the present and FEM 
solutions. 

Comparison of analytical solution with experimental 
results 

It is difficult to obtain accurate buckling loads in panels for 
buckling due to mechanical loading. They are mainly asso¬ 
ciated with obtaining the desired in-plane loading condi¬ 
tions (e.g., uniform stress), boundary conditions, and the 
flatness of the panels. In this study, all the equations were 
derived under the assumptions that the panel is perfectly 
flat, and the boundary edges are immovable and straight. In 
practice, wood-based panels have some imperfections, 11 and 
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boundary edges are movable because in-plane stresses act 
on the boundary edges. The hygrobuckling load or critical 
moisture content decreases as the degree of imperfection 
increases, 14,15 whereas it increases as the edge movement 
increases. In addition, the extent of the edge movement 
depends on in-plane stresses and panel geometries, includ¬ 
ing the core ratio. 1 Therefore, the presence of imperfections 
and the edge movement have combined effects on the buck¬ 
ling behavior when the panels are assembled in the core, 
which is one of the grid structures. It is not certain whether 
the boundary condition near surface panels is in a clamped 
or simply-supported condition even though surface panels 
of the hollow core door are fixed to the cores with adhesive. 
It is more likely in the clamped condition as Spalt and 
Sutton 11 assumed. 

By ignoring the presence of imperfections and edge 
movements, critical moisture contents calculated using Eq. 
(6) were 0.5% and 1.0% for hardboard and plywood, re¬ 
spectively. It was difficult to control changes in moisture 
content as much as is required, but the experimental results 
showed that hygrobuckling occurred when the change of 
moisture content was less than 2.0% for plywood and less 
than 1.0% for hardboard, as shown in Fig. 5b. 

The above results suggested that the present solu¬ 
tions are accurate for the all-edges-clamped condition. 
Therefore, the nonparameter hygrobuckling load and 
critical moisture content of plywood and hardboard were 
represented for further use as shown in Figs. 6 and 7, 
respectively. 3 

Equations (6) and (11) show that the critical moisture 
content depends on Poisson’s ratio, and it is inversely pro¬ 
portional to 1 + v for isotropic panels. It did not depend on 
the MOE at all for isotropic panels such as hardboard. 4 For 
orthotropic panels such as plywood, however, it depended 
on the MOE ratios. It should be noted in the case of ply¬ 
wood that the critical moisture content was much different 
for different directions even though the aspect ratio was the 


same (Table 6). To increase the critical moisture content at 
the same aspect ratio, therefore, the core should be placed 
perpendicular to the longitudinal direction of plywood. 


Conclusions 


Based on the displacement function and the energy method, 
buckling and hygrobuckling equations were derived for 
orthotropic panels subject to biaxial compression with all 
edges in both the clamped and simply supported conditions. 
The following conclusions were derived from this study. 

1. Using the displacement function with a double sine 
series, the results were similar to and simpler than 
those derived from Dickinson’s solution. 



Aspect ratio (a/b) 

Fig. 6. Nondimensional load parameter of orthotropic and isotropic 
panels as a function of aspect ratio. H = V 2I D U + ZD 66 



Fig. 5. Percent change of moisture content ( AMC ) and linear expan¬ 
sion ( LE ) under high relative humidity (RH) conditions (a) and center 
deflection with change of moisture content (b). a Solid triangles, hard¬ 
board MC; open triangles, plywood MC; solid circles, hardboard LE; 


open circles, plywood LE; solid down triangles, hardboard LE/MC; 
open down triangles, plywood LE/MC. b Solid line, hardboard - aver¬ 
age; dashed line, plywood-average; open circles, hardboard - observed; 
open squares, plywood - observed 
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Table 6. Comparison of critical moisture contents for different direc¬ 
tions of panels for the all-edges-clamped condition at the same aspect 
ratio 


Aspect ratio 

Plywood 

Hardboard 

alb 



(alh = 100) 

0.3 

1.5229 

0.5158 

3.0 

7.3658 

4.6321 

bla 



(, b/h = 100) 

0.3 

0.8146 

0.5056 

3.0 

13.7032 

4.6321 


Results (critical moisture contents) are percents 



Fig. 7. Changes of critical moisture content (MC„) by aspect ratio. 

Solid lines, plywood; dotted lines, hardboard; mangles, alh = 150; open 

circles, alh = 100; filled circles, alh = 50 

2. The derived equations could be applied to calculate the 
effects of material constants and the geometry of panels 
on hygrobuckling by ignoring the presence of imperfec¬ 
tions of the panel and the edge movement. 

3. The critical moisture content depended on Poisson’s 
ratio and was inversely proportional to 1 + v for iso¬ 
tropic panels. It did not depend on MOE at all for 
isotropic panels, whereas it did depend on the MOE 
ratios for orthotropic panels. 

4. The critical moisture content of plywood was twice as 
large as that for hardboard due to the difference in linear 
expansions between the two types of panel. 

5. Application of the optimum thickness and aspect ratios 
obtained by the derived equations could improve the 
hygrobuckling resistance without other chemical treat¬ 
ment to reduce the linear expansion of wood-based 
panels. It might be better to increase the aspect ratio 
rather than the thickness ratio (alh) from an economical 
viewpoint. For plywood, core ply should be placed 
perpendicular to the longitudinal direction. 


Appendix: abbreviations and symbols 



extensional and bending stiffness, respec¬ 
tively 

E, 

modulus of elasticity (MOE) in x and y 
directions, respectively 

G n 

modulus of rigidity in xy plane 

V- 

v l] 

Poisson’s ratio 

h 

panel thickness 

a, b 

panel length in x and y directions, respectively 

x, y, z 

coordinate system 

w 

displacement in z direction 

AMC 

moisture content change 

a„ a y 

coefficient of linear expansion in x and y 
directions, respectively 

LE 

linear expansion ( aAMC ) 

N„ N y 

resultant in-plane forces per unit length in x 
and y directions, respectively 

N?,N? 

hygroscopic forces in x and y directions, 
respectively 
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